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Abstract During surveys in central Florida of the zombie-ant fungus Ophiocordyceps camponoti-floridani, which 
manipulates the behavior of the carpenter ant Camponotus floridanus, two distinct fungal morphotypes were dis- 
covered associated with and purportedly parasitic on O. camponoti-floridani. Based on a combination of unique 
morphology, ecology and phylogenetic placement, we discovered that these morphotypes comprise two novel 
lineages of fungi. Here, we propose two new genera, Niveomyces and Torrubiellomyces, each including a single 
species within the families Cordycipitaceae and Ophiocordycipitaceae, respectively. We generated de novo draft 
genomes for both new species and performed morphological and multi-loci phylogenetic analyses. The macro- 
morphology and incidence of both new species, Niveomyces coronatus and Torrubiellomyces zombiae, suggest 
that these fungi are mycoparasites since their growth is observed exclusively on O. camponoti-floridani mycelium, 
stalks and ascomata, causing evident degradation of their fungal hosts. This work provides a starting point for 
more studies into fungal interactions between mycopathogens and entomopathogens, which have the potential to 
contribute towards efforts to battle the global rise of plant and animal mycoses. 
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INTRODUCTION 


Fungi occupy a wide array of ecological niches as decomposers, 
mutualists, and parasites of plants, animals and other fungi. 
Mycoparasites of other parasitic fungal lineages can impact 
ecosystem composition and disease dynamics by modulating 
their hosts' population size and transmission rate (Blackwell 
& Vega 2018). Despite these perceived ecosystem impacts 
and their biocontrol potential, mycoparasites are generally 
understudied, which is especially true for those that attack 
entomopathogenic fungi. 


The order Hypocreales contains the widest diversity of animal 
parasites among the kingdom Fungi. Most hypocrealean fungi 
are parasites of plants and arthropods, especially insects, 
although some species are known to parasitize spiders, 
nematodes, rotifers and even immunocompromised humans, 
as well as other fungi (Samson et al. 1988, Kepler et al. 2013, 
Lombard et al. 2015, Araujo & Hughes 2016). The genus Ophio- 
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cordyceps comprises approximately 300 species, strictly as- 
sociated with insects belonging to 13 orders (Crous et al. 2004, 
Araujo & Hughes 2016). Among these, one particular group 
stands out for its intriguing and bizarre biology, the Ophio- 
cordyceps unilateralis clade sensu Araüjo et al. (2018), which 
infect and manipulate the behavior of ants, mostly of the tribe 
Camponotini, across the globe (Andersen et al. 2009, Evans et 
al. 2011a, Araujo et al. 2015, 2018). Typically, Ophiocordyceps- 
infected ants, such as the Florida carpenter ant Camponotus 
floridanus, are manipulated to leave their nest and ascend 
vegetation, where they exhibit a fungus-adaptive 'death-grip' 
behaviour (Andersen et al. 2009, Araüjo & Hughes 2019, Will 
et al. 2020). Species within the O. unilateralis clade are highly 
specialised heterotrophs that are able to form epizootics, often 
infecting hundreds of ants within a small area of forest (Evans 
1982). After the spores encounter the host, penetrate and 
overcome its defences, the fungus proliferates as yeast-like 
cells in the haemocoel (see Araujo et al. 2020: f. 2d—f). Once 
established inside the host, the fungus produces secondary 
metabolites, proteases and other (small) secreted bioactive 
compounds to interact with its host and adaptively manipulate 
its behaviour (De Bekker et al. 2021). After the fungus kills the 
host, the yeast-like cells are converted into hyphae forming 
an endosclerotium, a compact mass of fungal mycelium that 
rapidly fills the host body after death (see Andersen et al. 2009: 
f. 3). Ophiocordyceps then utilizes the ant's body as a platform 
to grow the spore-producing structures needed for transmis- 
sion to the next host (Evans et al. 2011a, Hughes et al. 2011, 
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De Bekker et al. 2015, Araüjo et al. 2018). However, despite 
being sophisticated parasites themselves, but not unlike most 
(if not all) life on earth, Ophiocordyceps species are also 
parasitized by other hypocrealean fungi (Evans et al. 2011a, b, 
Andersen & Hughes 2012, Araujo et al. 2020). 


Mycoparasitism, or the parasitism of one fungus by another 
(Kirk et al. 2008), has independently and repeatedly appeared 
in a variety of fungal lineages along their evolution (Boddy 
2016, Herrera et al. 2016, Blackwell & Vega 2018). Thus far, 
mycoparasites associated with four other Ophiocordyceps spe- 


cies from China and Thailand have been described (Wang et 
al. 2015b, Zhong et al. 2016, Xiao et al. 2018). However, none 
ofthese species are associated with a behaviour-manipulating 
Ophiocordyceps. In addition, mycoparasites growing on other 
Ophiocordyceps-manipulated ants, such as Ophiocordyceps 
camponoti-rufipedis in Brazil (Evans et al. 2011a, b, Andersen 
& Hughes 2012) and Ophiocordyceps paltothyrei in Ghana 
(Araujo et al. 2020), have been reported. These records, along 
with our unpublished observations of mycoparasites on Ophio- 
cordyceps across North and South America, as well as in Africa 


granadensis on its host, Camponotus novogranadensis, in Atlantic rainforest, Itacolomi, Minas Gerais, Brazil (Note yellow perithecia on the subiculum in b; c—d. 
on O. camponoti-rufipedis on Camponotus rufipes in Atlantic rainforest, Viçosa, Minas Gerais, Brazil; e. torrubiellomyces-like perithecia on Ophiocordyceps 
camponoti-novogranadensis on, Camponotus novogranadensis, habitat as a, b; f. dark perithecia produced on the mycelium of Ophiocordyceps oecophyllae, 
on Oecophylla smaragdina in rainforest, Licuala State Forest, Queensland, Australia. 
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and Australia (JPM Araujo & HC Evans pers. obs., Fig. 1), sug- 
gest that the tri-trophic interactions that we report here are not 
unique to Florida. Instead, they are an example of a common 
worldwide phenomenon. Despite this, formal species descrip- 
tions and reports have remained limited: no genomes have 
currently been sequenced, and little research has been done on 
their biology and the effects that these mycoparasitic lineages 
have on Ophiocordyceps disease dynamics and transmission. 
Here, we describe two new genera, Niveomyces (Cordycipita- 
ceae) and Torrubiellomyces (Ophiocordycipitaceae) associated 
with and parasitic on Ophiocordyceps camponoti-floridani in 
Central Florida. Our proposal is supported by a polyphasic ap- 
proach combining morphological, ecological and phylogenetic 
data. While we predicted the placement of both morphotypes 
to reside in less data-rich parts of the hypocrealean tree (com- 
plicating culture identification through GenBank alignments) 
and obtaining PCRs from DNA extracted directly from field 
specimens proved to be difficult, we produced draft genomes 
of both species. These draft genomes were subsequently used 
to obtain sequences for multi-loci phylogenetic analyses and 
as alignment databases to identify the correct isolates obtained 
from additional specimens. These genomes also add consider- 
able data to the current low number of available mycoparasite 
genomes. 


MATERIALS AND METHODS 


Field sampling 


Ophiocordyceps camponoti-floridani-manipulated ant cadavers 
of C. floridanus with visible mycoparasitic growth of both mor- 
photypes were collected from the Black Hammock Wilderness 
Area (N28°42'04.7" W81°09'32.0") and Little Big Econ State 
Forest (N28°41'14.7" W81°09'33.4") in Central Florida. Two 
morphotypes were readily recognised in the field either by their 
characteristic cotton white hyphae that consistently covered 
and overgrew the host and Ophiocordyceps synnemata (Niveo- 
myces coronatus), or by the dark perithecia (Torrubiellomyces 
zombiae) that arose directly from the fungal host. Our collection 
permits were provided by the Seminole County's Leisure Ser- 
vices Department, Greenways and Natural Lands Division and 
the Florida Department of Agriculture and Consumer Service's 
Florida Forest Service. 


Fungal culturing 


To isolate N. coronatus, sterile water droplets of ~1 uL were 
pipetted onto parasitised synnema. Because of the high hydro- 
phobicity of the spore structures, droplets stayed intact and 
spores were released onto the water surface. Droplets with 
spores were streaked onto potato dextrose agar (PDA; BD 
Difco) and incubated at room temperature. After seven days, 
mycelium was transferred from single colonies to fresh PDA 
plates with a sterile inoculation loop for further isolation. To 
isolate T. zombiae, a sterile inoculation loop was used to pick 
up young, bright white, not-yet-matured fungal growth exhibited 
on top of O. camponoti-floridani and to inoculate PDA plates 
using the T-streak method. After incubating at room temperature 
for 14 d, single colonies were transferred to fresh PDA plates 
with a sterile inoculation loop for further isolation. Cultures of 
these isolates (ex-types) are deposited in the culture collec- 
tion of the Westerdijk Institute (CBS 149186 = BH-Nc-1D-3 for 
N. coronatus) and (CBS 149187 = BH-Tz-4E-4 for T. zombiae), 
respectively. 


To document culture characteristics of both species, the centre 
of fresh PDA plates was inoculated with mycelium from the two 
isolates using a sterile 4 mm-diam cork borer (Cole-Parmer). 
One half of the plates were incubated at room temperature 


and subject to daily light fluctuations in the lab. The other half 
was kept in the dark inside an incubator (Panasonic) kept at 
25 °C. After 6—8 wk of radial growth on PDA, the diameter of 
the cultures was measured and the mycelium was examined 
for presence of growth differentiation. 


To grow mycelium for DNA extractions, a flame-sterilised inocu- 
lation loop was used to scrape a small amount of mycelium from 
colonies growing on PDA plates to inoculate a sterile 250 mL 
Erlenmeyer flask containing 50 mL Sabouraud dextrose broth 
(SDB; BD Difco). Flasks were incubated at room temperature on 
a shaking platform (Fisher) at 120 rom. Mycelium was harvested 
from the Niveomyces liquid culture three days after inoculation 
while Torrubiellomyces was harvested after 9 d by pouring the 
culture over a Buchner funnel (Fisher) with Whatman Grade 
1 filter paper (Fisher) and applying suction. The mycelial- 
impregnated filter paper was pressed flat by hand between 
paper towels to remove any remaining liquid. A 1 cm? piece of 
the dried mycelium was placed in a 2 mL microcentrifuge tube 
(USA Scientific) containing two metal ball bearings (5/32" type 
2B, grade 300, Wheels Manufacturing) and snap-frozen in liquid 
nitrogen for tissue disruption and DNA extraction. 


Morphology 


To assess the macromorphological features, images were taken 
using a Canon EOS 7D Mark II camera fitted with a 35 mm 
lens. To investigate their micromorphological features, fungal 
tissues were mounted on microscope slides with a drop of either 
lactic acid, in the case of Torrubiellomyces, or lacto-fuchsin 
stain (0.1 g acid fuchsin powder and 100 mL 85 % lactic acid) 
for Niveomyces to aid visualisation of taxonomically informa- 
tive structures. The slides were visualized using a Leica DMi8 
inverted microscope, mounted with a Leica MC 170 HD camera 
(Leica Microsystems). Type materials (holotypes and paratypes) 
are deposited at the New York Botanical Garden Herbarium 
(type numbers NY4434800 and NY4434801 for N. coronatus 
and T. zombiae, respectively). 


DNA extraction, library preparation and whole-genome 
sequencing 


While the morphology of the fungi described in this study ap- 
peared to be unique compared to currently described species, 
we predicted that the sequence submissions of potentially relat- 
ed species would be vastly underrepresented in GenBank. This 
complicated direct genetic identification of sampled specimens 
and derived potential isolates based on PCR amplification. 
As such, draft genomes were generated for both new spe- 
cies using DNA extracted directly from collected specimens 
to obtain sequences for phylogenetic analyses and serve as 
genetic references to confirm the identity of isolated cultures 
through alignments. Under a dissecting microscope, tissues 
of the mycoparasites were removed while taking careful con- 
sideration not to include O. camponoti-floridani tissue. These 
tissues were surface sterilized in 70 % ethanol and placed into 
microcentrifuge tubes (USA Scientific) along with two metal ball 
bearings (5/32" type 2B, grade 300, Wheels Manufacturing), 
and snap-frozen in liquid nitrogen using a 1600 MiniG tissue 
homogenizer (SPEX) at 1300 RPM for 30 s to disrupt fungal cell 
walls. Genomic DNA was extracted using a previously described 
phenol-chloroform extraction protocol (Will et al. 2020), which 
was quantified with a Qubit Fluorometer (Thermo Fisher) and 
the Qubit dsDNA High Sensitivity Assay Kit (Thermo Fisher). 
Subsequently, DNA libraries were prepared with the Nextera 
DNA Flex Library Preparation Kit (Illumina) for sequencing on 
an Illumina MiSeq Sequencer to generate 2 x 250 bp paired- 
end reads with a 50x target coverage. 
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Simplicillium obclavatum CBS 311.74 


100 


77 


68 


87 


67 


————— 


100 U Simplicillium obclavatum CBS 311.74 

— — Simplicillium lanosoniveum CBS 704.86 
Simplicillium lamellicola CBS 116.25 
Simplicillium lamellicola CBS116.25 


Torrubiella wallacei CBS 101237 


1 


00 


Parengyodontium album CBS504.83 


Parengyodontium album CBS368.72 


— Parengyodontium album CBS504.83 


Flavocillium acerosum CBS 418.81 


89 Flavocillium primulinum JCM 18525 
Flavocillium primulinum JCM 18527 
[ei avo primulinum JCM 18526 


lavocillium sp. CBS 639.85 
e 


82 


71 


79 


67 


canicillium aranearum CBS 726.73a 
Hevansia cinereus NHJ 3510 
r——Hevansia novoguineensis NHJ 4314. 
levansia novoguineensis NHJ 13117 
Hevansia arachnophilus NHJ 10469 
94 Hevansia novoguineensis NHJ 11923 
vansia sp. BCC28584 
levansia nelumboides BCC41864 
Hevansia novoguineensis NHJ 13161 
Gibellula leiopus BCC 16025 
Gibellula sp. NHJ 13158 


7i Gibellula pulchra NHJ 10808 
[36 Gibellula sp. NHJ 5401 
71 Gibellula sp. NHJ 10788 
mog — Gibellula cf. alba NHJ 11679 
Gibellula longispora NHJ 12014 


Lecanicillium antillanum CBS 350.85 


Eng 


lyodontium rectidentatum CBS 641.74 


Engyodontium rectidentatum CBS 206.74 
Engyodontium parvisporum IHEM 22910 
Engyodontium aranearum 
Engyodontium aranearum CBS 309.85 


[—- Ascopolyporus vilosus ARSEF 6355 


—Ascopolyporus polychrous P.C. 546 

rL— Niveomyces coronatus Niveo 

100 Niveomyces coronatus NY04434800 
Niveomyces coronatus FieldW 


90 


99 


100 Pseudogibellula formicarum BCC81493 
Pseudogibellula formicarum BCC84257 
Pseudogibellula formicarum CBS871.72 
Pseudogibellula formicarum CBS433.73 


r——— Blackwellomyces minutus BCC88269 


Blackwellomyces calendulinus BCC68500 
Blackwellomyces calendulinus BCC68502 
1001, Blackwellomyces cardinalis OSC93610 


ob |" Blackwellomyces cardinalis OSC93609 


lackwellomyces aurantiacus BCC85060 
"H ig lackwellomyces aurantiacus BCC85061 
Blackwellomyces lateris MFLU18 0663 
Blackwellomyces roseostromatus BCC 91358 
Blackwellomyces roseostromatus BCC 91360 
Blackwellomyces roseostromatus BCC 91358 
Blackwellomyces roseostromatus BCC 91359 
Blackwellomyces pseudomilitaris BCC 2091 
Blackwellomyces pseudomilitaris BCC 73634 
+ Blackwellomyces pseudomilitaris TBRC 3662 
Blackwellomyces pseudomilitaris BCC 1919 
Lecanicillium fusisporum CBS 164.70 


'—Lecanicillium fusisporum CBS 164.70 


Lecanicillium psalliotae CBS 532.81 


66 <<" Lecanicillium psalliotae CBS 101270 


Liangia sinensis YFCC 3104 
Liangia sinensis YFCC 3103 
Hyperdermium piperis CBS 116719 


54 


—-Hyperdermium pulvinatum P.C. 602 


83 


5 


IE 


89 [— Akanthomyces attenuatus CBS 402.78 


Akanthomyces lecanii CBS 101247 
Akanthomyces coccidioperitheciata NHJ 6709 
r- Akanthomyces aculeatus HUA 186145 
84 1 — — —— Akanthomyces pistillariaeformis HUA 186131 
Akanthomyces tuberculatus OSC 111002 
Leptobacillium leptobactrum IRAN1230 
100 'Leptobacillium leptobactrum CBS 771.69 


31 


91 


— Leptobacillium muralicola 
Samsoniella Sp. spat 09-050 
[| Samsoniella sp. spat 09-051 
Samsoniella ramosa YFCC 6020 
jSamsoniella kunmingensis YHH 16002 
- Samsoniella cristata YFCC 6021 
83 Samsoniella tortricidae YFCC 6013 
Samsoniella antleroides YFCC 6113 
81 — Samsoniella lanmaoa YFCC 6148 
r Samsoniella alpina YFCC 5836 
92|| t Samsoniella cardinalis YFCC 5830 
Samsoniella hepiali YFCC 661 
Samsoniella yunnanensis YFCC 1527 
Beauveria malawaensis ARSEF 5413 
84|— Beauveria scarabaeicola ARSEF 5689 
r Beauveria caledonica ARSEF 2567 
6 Beauveria blattidicola MCA 1727 
- Beauveria pseudobassiana ARSEF 3405 
100r— — — Beauveira diapheromeriphila MCA 1557 
Beauveria neogryllotalpae MV2498 
— Beauveria locustiphila HUA 179218 
LT. Beauveria staphylinidicola ARSEF 5718 
Beauveria bassiana ARSEF 1564 
de Beauveria brongniartii BCC 16585 
Beauveria brongniarti ARSEF 617 
100; Cordyceps polyarthra MCA 996 
Cordyceps polyarthra MCA 1009 
Cordyceps exasperata MCA2155 
— Cordyceps exasperata MCA 2288 


100 


1e 


[co 


ER 


Qo. Cordyceps farinosa CBS 111113 
Cordyceps takaomontana MCA 1806 
10098 Cordyceps tenuipes OSC 111007 
Cordyceps sp. EFCC 2535 
r— Cordyceps coleopterarum CBS 110.73 
r— Cordyceps takaomontana BCC 12688 
88| || Cordyceps cicadae RCEF HP090724-31 
1 r- Cordyceps albocitrina Spat, 07-174 
48 Cordyceps bifusispora EFCC 5690 
—— Cordyceps ochraceostromata ARSEF 5691 
Cordyceps caloceroides MCA 2249 
97i Cordyceps oncoperae AFSEF 4358 
Cordyceps rosea spat. 09-053 
Cordyceps militaris OSC 93623 


92 


89| 


99.5 Cordyceps GE EFCC 5886 
100, —— —— Verticillium sp. CBS 102184 
51 — — ——— Cordyceps pruinosa ARSEF 5413 
r- Mariannaea pruinosa ARSEF 5413 
—— Cordyceps nimena oa EGS 38.166 
100 — Çordyceps fumosorosea CBS 337.52 
L— Cordyceps javanica CBS 134.22 


r— Cordyceps sp. RCEF HP090724-04C 
L- Cordyceps amoenerosea CBS 729.73 


Niveomyces 


Fig. 2 Maximum likelihood tree of Cordycipitaceae obtained with a concatenated dataset of SSU, LSU, TEF, RPB1 and RPB2. Niveomyces gen. nov. is indi- 
cated in bold font. The whole analysis tree of the order Hypocreales is depicted in the top-right corner, with the position of Cordycipitaceae highlighted in red. 
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Pleurocordyceps agarica YHHPA1303 


94 m Pleurocorc 


Pleurocorc 
Pleurocorc 


CT Pleurocordyceps agarica YHHPA1305 
Pleurocordyceps nipponica BCC18108 
Pleurocordyceps hipponica NHJ7727 
Pleurocordyceps nipponica BCC2325 
Pleurocordyceps nipponica NHJ4268 
Pleurocordyceps nipponica BCC1682 
Pleurocordyceps nipponica BCC1881 
Pleurocordyceps marginalradians MELU 17 1582 


fyceps marginaliradians MFLUCC 17 2276 


Pleurocordyceps phaothaiensis BCC84552 
Pleurocorüyceps phaothaiensis BCC84557 
Bleurocordyceps phaothaiensis BCCB4552 


ceps phaothaiensis BCC84553 


yceps phaothaiensis BCC84551 
Pleurocordyceps phaothaiensis BCCB4553 


Pleurocordyceps phaothaiensis BCC84551 
Pleurocordyceps aurantiaca MFLUCCTT 2114 
84 Bleurocordyceps aurantiaca MELUCCI7 1394 
Pleurocordyceps aurantiaca MFLUTI 1394 
Pleurocordyceps aurantiaca MFLU17 1393 
Pleurocordyceps aurantiaca MFLUCC17 2113 
Pleurocordyceps lianzhouensis GDGM 20918 
Pleurocordyceps lanzhouensis GIMMYOGO 
Pleurocordyceps ramosopulvinata EFCC5566 
Pleurocordyceps ramosopulvinata SU-65 
Pleurocordyceps sinensis CN80 2 
Pleurocordyceps sinensis ARSEF 1424 
Pleurocordyceps sinensis BL 4 
Pleurocordyceps sinensis GIMCC 3570 
Pleurocordyceps yunnanensis HHP 1006 
Pleurocordyceps yunnanensis YHCPY 1005 
Pleurocordyceps yunnaensis NBRC101760 
Pleurocordyceps yunnaensis NBRC 109984. 
Pleurocordyceps yunnaensis NBRC 109985 
Pleurocordyceps cf. yunnaensis EFCC 12075 
Polycephalomyces sp. NBRC 109987 
Polycephalomyces sp. NBRC110224 
Polycephalomyces sp. NBRC 109990 
Polycephalomyces sp. NBRC109988 
Torrubiellomyces zombiae NY04434801* 
Torrubiellomyces zombiae FieldB 
` Torrubiellomyces. zombiae Polyceph 
98 y Perennicordyceps prolifica INS-F 18547 
E Perennicordyceps pralifica TNS-F-18481 
Perennicordyceps cuboidea CEM 1514 


99 L Perennicordyceps cuboidea NBRC101740. 
96l 


82 


98 


86 


Torrubiellomyces) 


Perennicordyceps ryogamiensis 1017! 

Berennicordyceps ryogamiensis NBRC103837 

Perennicordyceps ryogamiensis NBRC103842 
Perennicordyceps paracuboidea NBRC101742 


Perennicon iboidea NBRC100942 
Cordyceps gi EE NERC AEN 
ENEE pleuricapitata NBRG100745 
Cordyceps pleuricapitata NBRČ100746, 

100 Polycephalomyces formosus CGMCC5.2208 
Polycephalomyces formosus CGMCCS 2006 
Polycephalomyces formosus CGMCCS 2204 
Polycephalomyces formosus CGMCC5.2207 
Polycephalomyces formosus CGMGC5.2205 


Polycephalomyces formosus CGMCC5.2203 
Harposporium harposporiferum ARSEF 5472 


53 Drechmeria zeosporum CBS 335.80 
83 86 [48 — — ——— Drechmeria sinensis CBS 567.95 
86 Drechmeria balanoides CBS 250.82 
86 Drechmeria sp. TNS F18495 
71 — — — — ~ Drechmeria gunnii OSC 76404 
Purpurèocillium atypicolum CBS 744.73 
Ke Purpureocillium átypicolum GNJO70713-08 
r Burpureocillium takamizusanense NHJ 3582 
43 Purpureocillium lilacinum CBS 284.36 
en = — Tolypocladiumm fractum OSC 110990 
ab Tolypociadium japonicum OSC 110991 
e {100 — lo ypocladium longisegmentum OSC 110992 
Tolypócladium capitatum OSC 71233 
9 Tolypócladium ophioglossoides CBS 100239 


rag — Tolypocladium cyli 
EK d 
Ophiocordyceps variabilis OSC 111003 

Ophiocordyceps formosana TNM F13893 
Ophiocorcyceps purpureostromata TNS F1843 


343 Ophiocordyceps konnoana EFCC 7 
phiocordyceps ravenelli OSC 110995 


bei L100- Ophiocordyceps nigrella EFCC 9247 
Ss Ophiocordyceps clavata CEM1762 
E ophiocorayceps entomorrhiza 16250 
Ophiacordyceps heteropoda EFCC 10125 


a 100 Ophiccordyceps sp. OSC 151905 
— phiocordyceps sp. OSC. 151904 


rosporur NRRL 28025 
OSC 71235 


77 


83 
-og — Ophiocordyceps amazonica HUA 186113 
9 Hoo E gracilis EFCC 8572 
Tog Ophiocordyceps gracilssima HUA 186132 
oaa 182 Ophiocort yceps melolonthae OSG 110993 
(TEE neovolkiana OSC 151903 
fgg Ophiocordyceps rhizoidea NHJ 12522 
Ophiocordyceps communis NHJ 12581 
57 Ophiocordyceps agriotidis ARSEF 5692 
Ophiocordyceps pruinosa NHJ 12994 
98 80 Ophiocordyceps cochlidiicala HMAS 199612 
Tt ophiocordyceps acicularis OSC 128580 


97, T Ophiocordyceps citrina TNS F1853 

56 pun Ophiocordyceps crinalis/acicularis TNS F T8550 
[54 81— — Ophiocordyceps stylophora OSC 111000 
[9^ — Ophiocordyceps sinensis EFCC 7287 


[gg ———— Ophiocordyceps elongata OSC 110989 
n 95 E Jeep? “Ophiocordyceps humbertii MF 16b 
6287 


E Ophiocordyceps tiputini QÓN 
s rr Cphiccoreyosps daceti MFO1 


is Ophiocordyceps monacidis MEA ` : 
68 r Ophiocordyceps ponerinarum HUA 186140 


L— Ophiocordyceps kniphofioides HUA 186148 
M — ——— ophiocordyceps oecophyllae OECO1 
36 — o Ophiocordyceps unilateralis s.l. Gh41 ; 
Ophiocordyceps rami MY6736 
Ophiocordycens camponoti-saundersi C40 


3 Ophiocordyceps ootakii J13 


32, Ophiocordyceps camponoti-leonardi C36 
BEE EE polyrhachis-furcata P39 


Ophiocordyceps halabalaensis MY 1308 
if ANT 


78 
Ophiocordyceops naomiplerceae DAWKS, 
100 73 E vede 


95 Ophiocordyceps blakebarnesii MISSOUS 
Ophiocordyceps pulvinata TNS-F 30044 
[————— Ophiocordyceps camponoti-balzani G104 d ti-floridani Flx1 (host 
a ft GEET EE canpgnoti-flondani E lxi (host) 
Las: Ophiocordyceps camponoti-rufipedis G108 
Ophiocordyceps camponoti-femorati FEMO2 
iocordyceps albacongiuae 
[pg — —  Ophlocordycens albacongl RC20, 2 
199 — —- Ophiocordyceps camponoti-nidulantis NIDUL2 
21 — — — — —— Ophiocordyceps camponoti reng jeri RENG2 
dz rog Ophiocordyceps Cam onoti-hippocrepidis HIP C 
3 Ophiocordyceps kimflemingiae SC30. 
Ophiocordyceps camponoti-bispinosi OBISS. 


Qphiocorayepes forquignonii OSC 151909 
[1 Gphiocordyceps forquignonii OSC 151908 
Ophiocordyceps nutans OSC 110994 
Ophiocordyceps lloydii OSC 151913 


49 
94 r Ophlocórd ene fulgoromorhla HUA 186142 
z EHS EES tulgoromorphila HUA 186139 
r=- Ophiocordyceps irangiensis OSC 128577 
100, Ophiocords ceps myrmecophila HMAS 199620 


63 90 


63 


94 14 


100 74 


1100 Ophiocordyceps irangiensis 128578 
Ophiocordyceps formicarum TNS F 18565 
700 ‘Ophiocordyceps sp.OSC 110998 
Ophiocordyceps sp. OSC 151909 
3ppCPhlocordyceps australis HUA 186147 
ši [oO Dphiocordyceps australis HUA 186097 
51 Ophiocordyceps buquetii HMAS 199613 
16 H SC 
Bd ‘Ophiocordyceps diabolica BDS32 
= ophiocordyceps palothyrel Patt 
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Fig. 3 Maximum likelihood tree of Ophiocordycipitaceae obtained with a concatenated dataset of SSU, LSU, TEF, RPB1 and RPB2. Torrubiellomyces gen. 
nov. is indicated in bold font. The host of Niveomyces and Torrubiellomyces, O. camponoti-floridani, is indicated in red. The whole analysis tree of the order 
Hypocreales is depicted in the top-right corner, with the position of Ophiocordycipitaceae highlighted in red. 
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Genome assembly and gene prediction 


Prior to assembling the two mycoparasite genomes, the raw 
sequence data were filtered and trimmed using the BBDuk 
plugin in Geneious Prime v. 20.2.3 with default parameters. 
Subsequently, to confirm the quality of these trimmed reads, 
fastQC was used (Andrews 2010). The genomes were then 
assembled de novo using the SPAdes assembly algorithm 
(Bankevich et al. 2012) and the quality of the assemblies was 
confirmed through QUAST (Gurevich et al. 2013). As expected 
from samples taken directly from the field, the QUAST outputs 
showed bacterial contamination in the genomic data, indicated 
by the presence of two distinct mean G-C 96 peaks; one large 
peak comprising fungal reads (30—80 96) and a second, much 
smaller peak comprising bacterial contaminants (0—30 96). The 
bacterial contaminants were removed from both genomes by 
manually removing sequences with a mean G-C % that fell 
within the contaminant peak (i.e. « 30 96). The effectiveness of 
this bacterial filtering was confirmed using MG-RAST (Keegan 
et al. 2016). We determined the completeness of the genomes 
after bacterial filtering with BUSCO (Seppey et al. 2019), using 
the Hypocreales lineage (fungi odb9), and CEGMA (Parra et al. 
2007). Scaffolds shorter than 1000 bp were discarded. Gene 
predictions were performed with Augustus v. 3.3 (Stanke & 
Morgenstern 2005) using the previously generated parameters 
for O. camponoti-floridani (Will et al. 2020) and the software 
parameters for Fusarium graminearum. The draft genomes 
are available through GenBank under the accession num- 
bers: JADHZA000000000 (Torrubiellomyces zombiae) and 
JAFEME000000000 (Niveomyces coronatus). The genome 
assembly, gene predictions and functional annotations can 
also be interactively analysed and downloaded through https:// 
fungalgenomics.science.uu.nl. 


Functional annotation 


The predicted proteins in our draft genomes were functionally 
annotated using PFAM (El-Gebali et al. 2019) and mapped to 
their corresponding gene ontology (GO) terms. Transmembrane 
domains were annotated using TMHMM v. 2.0 (Krogh et al. 
2001) and signal peptides using SignalP-5.0 (Nielsen et al. 
1997). Proteins with a secretory signal were considered small 
secreted proteins if they were shorter than 300 amino acids and 
did not contain a transmembrane domain (except in the first 
40 amino acids). We used a SMURF-based pipeline to predict 
secondary metabolite clusters (Khaldi et al. 2010, De Bekker 
et al. 2015). BlastP, with an E-value cutoff of 1e-10, was used 
to search the MEROPS database for proteases (Rawlings et 
al. 2018). 


DNA extraction, PCR and phylogenetics 


Using our de novo assembled and annotated mycoparasite 
genomes as protein Blast databases, we obtained and verified 
sequences for ribosomal 18S (SSU), ribosomal 28S (LSU), 
translation elongation factor 1-alpha (TEF), and RNA Polymer- 
ase Il Subunits (RPB1 and RPB2) for phylogenetic placement. 
We also extracted DNA from additional specimens and liquid 
cultures of the isolates for PCR amplification of these genes. 
We extracted DNA as previously described (Will et al. 2020). We 
amplified genes with Phusion polymerase (New England Bio- 
labs) using the primers and PCR programs published in Araujo 
et al. (2018), with cycle lengths and temperatures adjusted as 
per the recommendations provided in the Phusion polymerase 
protocol. We aligned the obtained sequences in Geneious to 
a database comprised of 531 species (Table 1) that broadly 
represented the order Hypocreales. Each locus was individually 
aligned with MAFFT (Katoh & Standley 2013) and concatenated 
into a single combined dataset using Geneious v. 11.1.5. The 
concatenated files, along with a position (POS) file for each 


gene, were imported into CIPRES (Miller et al. 2012). The final 
alignment length was 4770 bp: 1244 bp for SSU, 939 bp for 
LSU, 963 bp for TEF, 639 bp for RPB1 and 985 bp for RPB2. 
We performed maximum likelihood analysis with RAXML v. 8.2.4 
(Stamatakis 2014) on a concatenated dataset containing all five 
genes. The dataset consisted of 11 partitions, two for SSU and 
LSU and nine for each codon position of the three protein coding 
genes: TEF, RPB1 and RPB2. We employed the GTRGAMMA 
model of nucleotide substitution during the generation of 1000 
bootstrap replicates. Visualization and graphic adjustments 
were made in Dendroscope (Huson & Scornavaca 2012) and 
further edited in Adobe Illustrator. 


RESULTS 


Phylogenetics 


The phylogenetic results recovered the overall topology pre- 
sented in previous studies (Quandt et al. 2015, Kepler et al. 
2017, Araujo et al. 2018, Araujo & Hughes 2019). To determine 
the phylogenetic placement of the two mycoparasitic species, 
a comprehensive phylogenetic tree of the order Hypocreales 
was generated, adapted from Araujo & Hughes (2019). Based 
on the phylogenetic results, both species, N. coronatus and 
T. zombiae, formed distinct and well-supported monophyletic 
clades, BS = 100 for Niveomyces and BS = 98 for Torrubiel- 
lomyces (Fig. 2, 3). 

According to the data, Niveomyces occupies a basal branch 
within the family Cordycipitaceae (Fig. 2), while Torrubiel- 
lomyces sits in a basal clade within Ophiocordycipitaceae 
(Fig. 3). Niveomyces formed a unique, distinctive and relatively 
long-branched clade, while Torrubiellomyces fell within the 
Polycephalomyces clade. In order to investigate the relation- 
ships of Torrubiellomyces, we sampled Polycephalomyces 
s.lat. (Polycephalomyces insertae sedis, Pleurocordyceps and 
Perennicordyceps); including a range of species representing 
distinct ecologies, such as animal and fungal parasites (Kepler 
et al. 2013, Matoéec et al. 2014, Xiao et al. 2018, Wang et al. 
2021) (Fig. 3). Our phylogeny suggests the Polycephalomyces 
clade as the most basal lineage within Ophiocordycipitaceae 
(BS = 53) and strongly supports Torrubiellomyces as a distinct 
genus (BS = 98), closely related to a clade strictly associated 
with insects: Perennicordyceps cuboidea, Pe. Paracuboidea, 
Pe. ryogamiensis (all on Coleoptera), and Pe. prolificus (on 
Hemiptera) (Fig. 3) (Kepler et al. 2013). 


Taxonomy 


Based on a combination of morphological, ecological and phylo- 
genetic data, we introduce two new genera and two new spe- 
cies of mycoparasites within the Hypocreales. Torrubiellomy- 
ces zombiae and Niveomyces coronatus were both collected 
parasitizing Ophiocordyceps camponoti-floridani, a ubiquitous 
entomopathogen of the ant Camponotus floridanus in Florida, 
USA. 


Niveomyces J.P.M. Araüjo & C. de Bekker, gen. nov. — Myco- 
Bank MB 839229 


Etymology. Name reflects the 'snowy' (Lat.: niveus) appearance of this 
fungus. 


Type species. Niveomyces coronatus J.P.M. Araüjo & C. de Bekker 


Diagnosis: Niveomyces is diagnosed by its mycoparasitic nature, the produc- 
tion of spiky, white, slender, velvety synnemata and unique characters of 
the conidiogenous cells, which exhibit multiple denticles along the phialides 
with a crown-like apex, producing conidia singly. 


Mycelium white to pale yellow, often covering the host entirely. 
Vegetative hyphae septate and hyaline. Synnemata multiple, 
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Fig. 4 Niveomyces coronatus growing on Ophiocordyceps camponoti-floridani, a pathogen of the ant Camponotus floridanus. a. View of the tri-trophic system 
ant-entomopathogenic fungi-mycoparasite; b. close-up of N. coronatus synnemata; c. PDA culture after 60 d; d. close-up of culture edge; e. close-up of sporo- 
dochia formed in culture; f. layer of phialides (hymenium); g. close-up of apical and lateral conidiogenous cells; h. conidium. — Scale bars: f-g = 10 um, h = 5 um. 
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Spiky, erect, slightly sinuous to straight, not branched, tapering 
towards the apex, covered by hymenium-like layer of conidio- 
genous cells. Conidiogenous cells polyblastic, elongated, ir- 
regular, hyaline, cylindrical, with characteristic denticles that are 
crowded on the apical part and less frequent towards the 
base. Conidia globose to ovoid formed singly on the denticles. 
Produces micromorphological features, such as conidiogenous 
cells and conidia, identical on both specimen and in culture. 
Sexual morph unknown. 

Hosts — Entomopathogenic fungi. 

Distribution — USA, but probably worldwide. 


Niveomyces coronatus J.P.M. Araüjo & C. de Bekker, sp. nov. 
— MycoBank MB 844049; Fig. 4 


Etymology. Name reflects the characteristic crown of denticles on top of 
the conidiogenous cells. 


Typus. USA, Florida, Seminole County, Oviedo, Little Big Econ State For- 
est, N28°41'14.7" W81°09'33.4", over-growing Ophiocordyceps camponoti- 
floridani, a fungal pathogen of Camponotus floridanus, 10 June 2017, de 
Bekker (holotype NY04434800). 


Diagnosis: White mycelium covering the host almost entirely, producing 
multiple spike-like synnemata; exhibiting abundant characteristic coni- 
diogenous cells bearing multiple denticles, especially at the apical part 
producing globose to ovoid conidia. 


Mycelium white to light yellow, growing abundantly on the host. 
Synnemata multiple, not branching, white, slender, erect, arising 
from the subiculum that covers the host almost entirely, nar- 
rowing towards the end, averaging 311.9 x 65.2 um; covered 
by ahymenium-like layer of dense conidiogenous cells. Conidio- 
genous cells (12—)17(—25) x 1.5—2 um, polyblastic, cylindrical, 
hyaline, irregular, sometimes capitate, bearing crowded hyaline 
denticles on the apical part, often descending sparsely along al- 
most the entire cell. Conidia (3.7—)4.5 (-5.5) x 1.5-2(-2.5) um 
formed singly on the denticles, solitary, ovoid to globose, one- 
celled, hyaline and smooth-walled. 

Culture characteristics — Colonies on PDA reach a dia- 
meter of 65 —70 mm after 6 wk incubation at room temperature. 
Mycelium white during early stages becoming light yellow with 
age; remaining thin without spore formation when grown at 
25 °C in total darkness; at room temperature, under regular 
day-night light fluctuations, aerial hyphae and conidia similar 
to those observed on field specimens formed at the periphery 
of the colony after 8 wk. Synnemata were produced after 10 wk 
of incubation. 

Host — Ophiocordyceps camponoti-floridani. 

Distribution — Florida, USA. 

Ecology — Parasitic on Ophiocordyceps camponoti-floridani, 
an entomopathogen of the Florida carpenter ant, Camponotus 
floridanus, a ground-dwelling ant, commonly nesting in dead 
wood or soil. Ophiocordyceps-manipulated C. floridanus are 
predominantly found in elevated position, ranging from 0.1 m 
up to at least 2.5 m in height, clinging on and biting on epi- 
phytic plants of the family Bromeliaceae (specifically, Tillandsia 
recurvata and T. usneoides) in mesic hammock habitats with 
evergreen canopy. 


A similar species has been found parasitising both Ophio- 
cordyceps camponoti-novogranadensis on its host, Campo- 
notus novogranadensis, and O. camponoti-rufipedis on Cam- 
ponotus rufipes in remnant Atlantic rainforest in Minas Gerais, 
Brazil (H.C. Evans pers. obs.). Here, they are common myco- 
parasites — especially in the wet season — and may be exert- 
ing some control of Ophiocordyceps infections (see Fig. 1). It 
remains to be determined if these are N. coronatus or another 
member of the genus. 


Torrubiellomyces J.P.M. Araujo & C. de Bekker, gen. nov. — 
MycoBank MB 844048 


Etymology. Name reflects the resemblance to torrubiella-like fungi with 
the production of superficial perithecia on a subiculum directly on the host, 
without a stipe. 


Type species. Torrubiellomyces zombiae J.P.M. Araüjo & C. de Bekker 


Diagnosis: Torrubiellomyces is diagnosed by its mycoparasitic habit, its ini- 
tially white, then brown to black perithecia produced directly on the fungal 
host tissue. In culture, it exhibits a characteristic production of viscous 
conidial masses, similar to those produced by Pleurocordyceps, but without 
a stipe (synnema) on the host or in culture (Fig. 5, 6). 


Mycelium initially white, then brown growing on the insect 
sutures where the fungal host emerges; producing solitary or 
clusters of superficial perithecia that turn black and powdery 
with age; no synnemata or stalk produced on the host or in 
culture. Ascospores filiform, disarticulating into part-spores 
upon maturity. In culture, producing light cream mycelium with 
pools or masses of viscous conidia in sporodochia (Fig. 6c), 
never forming synnemata. Only one type of conidiogenous 
cell observed. 

Hosts — Entomopathogenic fungi. 

Distribution — USA, but probably worldwide. 


Torrubiellomyces zombiae J.P.M. Araüjo & C. de Bekker, sp. 
nov. — MycoBank MB844050; Fig. 5, 6 


Etymology. Name reflects its pathogenic association with the zombie-ant 
fungus O. camponoti-floridani, which manipulates its host's behaviour during 
infection. 


Typus. USA, Florida, Seminole County, Oviedo, Little Big Econ State 
Forest, N28°41'14.7" W81°09'33.4", on Ophiocordyceps camponoti-floridani 
pathogenic on Camponotus floridanus, 12 June 2017, de Bekker (holotype 
NY04434801). 


Diagnosis: Multiple perithecia produced solitarily or forming aggregations on 
multiple parts of the host mycelium, turning black and powdery with age. 
Asexual morph only observed in culture, forming pools of viscous conidial 
masses in sporodochia, lacking a stipe, one type of conidiogenous cell 
observed; producing ovoid, hyaline, smooth conidia. 


Mycelium scarce, pale to brown; growing sparsely, producing 
brown perithecia becoming black and powdery with age, pro- 
duced directly on the host synnemata or on host mycelium that 
emerges from the ant joints and sutures. Perithecia superficial, 
ovoid, solitary or forming dense aggregations, (420—)480(—525) x 
(205—)280(—310) um, rugose when dry. Asci hyaline, capi- 
tate, cylindrical, 225 x 6.5-7 um. Ascospores arranged in a 
spiral within the ascus, readily breaking into part-spores often 
still within the ascus. Part-spores hyaline, cuboid to globose, 
and often varying in size within a single ascus when mature, 
exhibiting a corn cob-like aspect (Fig. 5f), (1.8—)2.4(-3.2) x 
(1.5—)2(-2.6) um. Asexual morph not observed on any of the 
field-collected specimens. 

Culture characteristics — Colonies on PDA reaching a dia- 
meter of c. 40—48 mm after 6 wk incubation at room tem- 
perature. Mycelium white during early stages; becoming light 
cream with age, dense and reverse brown. Synnemata never 
observed, but conidial masses produced after 6—8 wk directly 
on the subiculum; pale cream, usually surrounded by sterile 
perithecium-like structures. Phialides producing large number of 
conidia, forming viscous masses (Fig. 6c); cylindrical to sub- 
ulate, usually slightly curved, (5.5—)8(—12) x (1.5—)1.7(-2.2) 
um, tapering gradually towards the apex. Conidia solitary, single 
celled, smooth-walled, hyaline, ovoid, 2.5—3 x 1.5—1.8 um. 

Host — Ophiocordyceps camponoti-floridani. 

Distribution — Florida, USA. 

Ecology — As above. 
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Fig. 5 Torrubiellomyces zombiae growing on Ophiocordyceps camponoti-floridani, a pathogen of the ant Camponotus floridanus. a. Close-up showing 
perithecia emerging from the antennal plate of the ant; b. general overview of a typical perithecial arrangement in clusters or less often singly produced; c—e. 
perithecia; f. mature ascospores disarticulated prior to release, forming corn cob-like asci; g. cluster of part-spores; h. ascus showing the ascospores already 
sub-divided into part-spores. — Scale bars: c—d = 100 um, e = 50 um, f = 10 um, g-h = 5 um. 


A similar species has been found parasitising O. camponoti- 
novogranadensis on its host, Camponotus novogranadensis 
in Brazil and O. oecophyllae infecting Oecophylla smaragdina 
in Australis (see above; Fig. 1e and 1f, respectively). Similar 
perithecia have also been found on the mycelium of Ophio- 
cordyceps oecophyllae (Araüjo et al. 2018), a pathogen of the 
weaver ant Oecophylla smaragdina in the rainforest of tropical 
Queensland, Australia. It remains to be determined if these 
records are T. zombiae or a related species of the genus. 


Draft genomes 


As part of the species descriptions, the genomes of the two 
novel fungal species, N. coronatus and T. zombiae, were se- 
quenced and assembled de novo. Our assembly of the N. coro- 
natus genome resulted in a total genome size of 31.95 Mbp, 
made up by 1357 contigs, with an N50 of 49324 bp and a 
G+C content of 51.11 %. The genome contains 8 930 protein 
encoding sequences of which 6433 sequences (72.04 96) 
were functionally annotated with PFAM domains and 4 159 
(46.57 %) received GO annotations. In addition, this genome 
is predicted to contain 766 genes with secretion signals, 1714 
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Fig. 6 Torrubiellomyces zombiae in culture. a. PDA plate after 60 d, arrows indicate the pools of viscous conidia produced in sporodochia scattered over the 
plate; b. overview of the specimen, white arrows indicate T. zombiae perithecia emerging from the fungal host tissue; c. close-up of sporodochium with pools 
of viscous conidia; d. conidia; e. cluster of phialides; f—I. phialides. — Scale bars: c = 1000 um, d = 3 um, e = 20 um, f-17 5 um. 


genes with transmembrane domains, 304 transcription factors, 
336 proteases and 53 secondary metabolite clusters. The set of 
gene predictions was determined to be 93.1 % complete (i.e., 
BUSCO completeness, CEGMA was 97.16 96). The assembled 
T. zombiae genome was calculated to be 27.11 Mbp in size, 
consisting of 1725 contigs with an N50 of 41529 bp and a G+C 
content of 51.49 %. Our annotation predicted the genome to 
contain 8 422 protein encoding sequences with 5 985 (71.06 96) 


of those sequences containing known PFAM domains and 
3812 (45.26 95) receiving GO annotations. Additionally, this 
genome was functionally annotated to have 675 genes with 
secretion signals, 1496 genes with transmembrane domains, 
242 transcription factors, 274 proteases and 32 secondary 
metabolite clusters. The set of gene predictions was found to 
be 94.48 96 complete (i.e., BUSCO completeness, CEGMA 
was 98.03 96) (Table 2). 
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Table 2 Draft genome details and statistics on genome assembly, gene predictions, annotations and completeness. 


Property 


Torrubiellomyces zombiae 


Niveomyces coronatus 


Sequences in assembly 

Total assembly length (Mbp) 
Assembly GC content (96) 

Assembly gaps (96) 

L50 number (£) 

N50 length (bp) 

Genes 

Gene length (median) 

Transcript length (median) 

Exon length (median) 

CDS length (median) 

Protein length (median) 

Spliced genes (total, 96) 

Exons per gene (median) 

Intron length (median) 

Introns per spliced gene (median) 
Gene density (genes / Mbp) 

Coding content of assembly (bp, 96) 
Proteins with internal stops (total, 96) 
Unique PFAM domains 

Genes with PFAM (total, 96) 

Genes with GO (total, 96) 

Genes with signalP (total, 96) 

Genes with TMHMM (total, 96) 
Genes annotated as TF (total, 96) 
Genes annotated as MEROPS protease (total, 96) 
Genes annotated as CAZyme (total, 96) 
Secondary metabolite clusters 
CEGMA completeness (96) 


1725 
27,11 

51,49 

0 

189 

41529 

8422 

1382 

1233 

294 

1230 

410 

6294 (74.73%) 
2 

64 

2 

310,68 
12740001 (47.0%) 
0 (0.0%) 

3803 

5985 (71.06%) 
3812 (45.26%) 
675 (8.01%) 
1496 (17.76%) 
242 (2.87%) 
266 (3.16%) 
224 (2.66%) 
32 

98,03 


1357 
31,95 

51,11 

0 

200 

49324 

8930 

1452 

1302 

323 

1299 

433 

6450 (72.23%) 
2 

67 

2 

279,49 
13987488 (43.78%) 
0 (0.0%) 

3786 

6433 (72.04%) 
4159 (46.57%) 
766 (8.58%) 
1714 (19.19%) 
304 (3.4%) 
319 (3.57%) 
296 (3.31%) 
38 

97,16 


BUSCO2 completeness (fungi_odb9) 


Complete: 94.48 % (Single-copy: 93.79 %, 
Duplicated: 0.69 %), Fragmented: 3.79 %, 
Missing: 1.72 % 


Complete: 93.1 % (Single-copy: 92.07 %, 
Duplicated: 1.03 %), Fragmented: 2.41 %, 
Missing: 4.48 % 


DISCUSSION 


Few mycoparasitic species of entomopathogenic fungi have 
been formally described, thus far, including recent records 
of Polycephalomyces on Ophiocordyceps species from Asia 
(Wang et al. 2015b, Zhong et al. 2016, Xiao et al. 2018). How- 
ever, mycoparasites associated with behaviour-manipulating 
Ophiocordyceps have only been noted in the field as an 
ecological aspect of those interactions (Andersen & Hughes 
2012, Araüjo et al. 2020, Mongkolsamrit et al. 2021). In this 
study, we present two new genera Niveomyces and Torru- 
biellomyces, which were recorded consistently infecting the 
zombie-ant fungus O. camponoti-floridani. We also provide 
their annotated draft genomes, which we used as a means 
to genetically identify the correctly cultured isolates that were 
obtained from field specimens and to extract sequences for 
phylogenetic analysis. Only a limited number of hypocrealean 
mycoparasites has been sequenced so far, including Tolypocla- 
dium ophioglossoides (Quandt et al. 2015), several Trichoderma 
species such as Trichoderma virens (Kubicek et al. 2011), 
Trichoderma atroviride (Kubicek et al. 2011) and Trichoderma 
reesei (Martinez et al. 2008), Escovopsis weberi (De Man et al. 
2016), and Clonostachys rosea (Karlsson et al. 2015). The draft 
genomes that we generated for N. coronatus and T. zombiae, 
therefore, represent a significant contribution to the still scarce, 
existing mycoparasite genomics data. Currently, these data are 
too scattered across the Hypocreales to conduct meaningful 
comparative genomics analyses into mycoparasite signatures. 
However, we hope that this study will stimulate more research 
into mycoparasitism and generate additional draft genomes in 
order to make such analyses more worthwhile. 


Both mycoparasites are well-supported by the comprehensive 
phylogeny as novel and unique lineages in the families Cordy- 


cipitaceae and Ophiocordycipitaceae of the order Hypocreales. 
Torrubiellomyces zombiae is placed as a new lineage within 
the Polycephalomyces clade, as sister to Perennicordyceps, 
a genus composed exclusively of entomopathogenic species. 
This suggests that its origins are from an insect-associated 
ancestor shared with Perennicordyceps. However, other spe- 
cies within the Polycephalomyces clade are also parasitic on 
entomopathogenic fungi. For example: Pleurocordyceps yun- 
naensis on Ophiocordyceps nutans, a pathogen of stink bugs; 
PI. aurantiaca on O. barnesii, a pathogen of melolonthid larvae 
and PI. agarica on an unidentified Ophiocordyceps species also 
pathogenic on melolonthid larvae; demonstrating the affinity of 
this group to exploit entomopathogens (Wang et al. 2015b, Xiao 
et al. 2018). Furthermore, Niveomyces coronatus resides within 
a part of the Hypocreales tree that, thus far, largely contains 
entomopathogens. This suggests that its mycoparasitism may 
have evolved from a previous animal parasitic relationship. 


Regarding morphological features, Niveomyces coronatus 
exhibits snow-white mycelium that often completely covers the 
host, producing multiple synnemata on a subiculum or directly 
on the host tissue; while its sister genus — Pseudogibellula — 
has gibellula-like conidiophores but which, unlike the phialidic 
heads of Gibellula, has heads of conidiogenous cells producing 
conidia sympodially on minute denticles, leaving protuberant 
scars (Samson & Evans 1973, Araüjo et al. 2020, Mongkolsam- 
rit et al. 2021). Pseudogibellula is a monotypic genus and the 
type species, P. formicarum has been described as a 'strongly 
competitive fungus on insect substrates and frequently exploits 
ant cadavers killed by other fungal pathogens' (Samson & 
Evans 1973). However, it was also reported to cause local 
epizootics on at least six ant species in evergreen forest in 
Ghana, as well as being a pathogen of several Homopteran 
hosts in cocoa farms (Samson & Evans 1973). Thus, P. for- 
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micarum has been considered — seemingly, ambiguously — as 
both an opportunistic mycoparasite and as an entomopathogen. 
Both Araujo et al. (2020) and Mongkolsamrit et al. (2021) have 
recorded P. formicarum in a purported mycoparasitic associa- 
tion with Ophiocordyceps pathogens of ants and leafhoppers, 
respectively. However, these authors reasoned that the fungus 
may also be an entomopathogen, based on evidence of primary 
infection of insect hosts (Homoptera: Cicadellidae) in both 
Brazil and the USA. In the latter, this involved both in vivo and 
in vitro studies of the interaction and from the results of patho- 
genicity experiments, it was concluded that P. formicarum is a 
primary entomopathogen and responsible for field epizootics 
ofthe glassy-winged sharpshooter, Homalodisca coagulata: a 
major agricultural pest in Florida (Kanga et al. 2004, Boucias 
et al. 2006). In fact, Boucias et al. (2006) noted differences in 
conidiophore morphology in the Florida isolate and considered 
that this could be a novel species of Pseudogibellula. A more 
detailed molecular analysis of the various geographical and 
host isolates of P. formicarum seems warranted, especially to 
compare the purported mycoparasitic and entomopathogenic 
strains. 

This apparent inter-kingdom jump — with mycoparasitic spe- 
cies evolving from entomopathogens to parasitise related 
entomopathogenic fungi — is analogous to mycoparasitism in 
the Urediniomycetes, where rust relatives of the genus Tuber- 
cularia are parasitic on rust fungi; having evolved from a plant 
parasitic lineage (Lutz et al. 2004a-c). The plant parasitic 
Helicobasidium sexual morph has a wide host range, whereas 
the Tubercularia mycoparasitic species show a high degree of 
specificity within their rust hosts. It remains to be confirmed if the 
mycoparasites described here, as well as P. formicarum s.lat., 
have similar levels of specificity within their entomopathogenic 
fungal hosts. 


Torrubiellomyces is only known from its sexual morph and it is 
easily recognized in the field by the formation of single or clus- 
ters of brown to black superficial perithecia that are produced 
directly on the mycelium of the Ophiocordyceps host fungus. In 
culture, Torrubiellomyces forms viscous conidia that are charac- 
teristic of species belonging to the Polycephalomyces clade. 
However, it differs from other closely related genera (Peren- 
nicordyceps, Pleurocordyceps and Polycephalomyces s.str.) 
by the absence of a stipe (synnema) supporting the formation 
of viscous conidia, which in Torrubiellomyces are produced in 
sporodochia in vitro (see Fig. 6c). 


Field observations 


The fungi that we consistently find growing on Ophiocordyceps 
have always been considered as their associated mycopara- 
sites. However, one could perhaps argue that they could also 
be growing saprophytically on dead insect or fungal tissue 
or act as entomopathogens that co-infect the insect host of 
Ophiocordyceps. The latter is especially enticing considering 
the phylogenetic placement of both species among fungal 
groups that broadly include entomopathogenic species. 
However, based on their unique morphology, which indicates 
that both N. coronatus and T. zombiae only grow on top of 
Ophiocordyceps tissue, it is more likely that they are indeed 
mycoparasites as posited in other studies dealing with interac- 
tions between entomopathogenic fungi and other antagonistic 
fungi (Wang et al. 2015b). Our field observations supports 
this conclusion. Although in some cases N. coronatus and 
T. zombiae growth was observed less than a week after a new 
O. camponoti-floridani-infected ant cadaver was found, they 
also appeared up to nine months after initial Ophiocordyceps 
infection, suggesting that these mycoparasites are able to in- 
fect Ophiocordyceps species at any stage of its development. 
Moreover, O. camponoti-floridani mycelium begins to emerge 


from a fresh ant cadaver within one to two days after its death, 
providing already sufficient tissue for the mycoparasite to in- 
habit. This seemingly rules out the possibility of either species 
being an entomopathogen that coinfects the ant. Moreover, 
already within the first hours after an ant's manipulation and 
following death, Ophiocordyceps completely colonizes the ant's 
body, consuming all ant tissue, besides the cuticle, to gain 
the energy needed to grow the fruiting body (De Bekker et al. 
2015). This makes it unlikely that either of the two species that 
we describe here would be able to saprophytically consume the 
insect cadaver. Moreover, the appearance of both mycoparasitic 
species in the first few weeks after death of the ant host, when 
O. camponoti-floridani is fresh and actively growing, seems to 
rule out the possibility of either species being an opportunistic 
saprophyte that merely feeds on dead O. camponoti-floridani 
tissue. In fact, even months after ant manipulation and death 
have taken place, Ophiocordyceps is often found to be alive 
andis still able to produce stalks with ascoma and consequently 
release spores. Taken together, our morphological and field 
observations confirm previous assumptions that the fungi we 
found in association with Ophiocordyceps species are indeed 
mycoparasites. 


Conclusions 


Here, we describe two new genera of fungi — Niveomyces and 
Torrubiellomyces — parasitic on the zombie-ant fungus, Ophio- 
cordyceps camponoti-floridani, from a small sample area in 
central Florida. Collections of similar fungi have been made 
from other Ophiocordyceps species, especially those attack- 
ing ants, in South America, Africa, Asia and Australia. It is 
likely, therefore, that such mycoparasites are pantropical and 
that these tri-trophic interactions are important contributors 
to the 'natural balance' in their respective ecosystems. The 
diversity and host specificity within these new genera, as well 
as in related genera, such as Pseudogibellula, remains to be 
determined. However, it is possible that such mycoparasitic 
genera harbour a potentially large and untapped reservoir of 
undocumented fungal diversity, especially when we consider 
the diversity of entomopathogenic fungi worldwide (Araujo A 
Hughes 2016). This is part of the 'hidden fungal biodiversity' 
described by Blackwell & Vega (2018) and evidence suggests 
that there was a diverse range of mycoparasites in existence 
over 400 million years ago (Berbee et al. 2017, Krings et 
al. 2017). Mycoparasites of plant pathogens also constitute 
this 'cryptic' biodiversity and are now being more intensively 
studied because of their potential for biological control of plant 
diseases. The classical biological control approach, involving 
surveys in the centres of origin or diversity of the target plant 
pathogen, has yielded a surprising diversity from relatively 
small sample sizes. Mycoparasites associated with frosty pod 
of cacao (Moniliophthora roreri) on its wild Theobroma host in 
the forests of western Ecuador (Evans et al. 2003) and those 
associated with coffee leaf rusts (Hemileia spp.) on their wild 
Coffea hosts in Africa (Colmán et al. 2021, Rodríguez et al. 
2021) provide evident examples. Recent estimates of extant 
fungal species, based on a fungal census of soils and a fungal/ 
plant ratio of 17/1, has put the number near six million spe- 
cies (Taylor et al. 2014). These authors concluded that: “98 % 
of fungi remain undescribed and that many of these species 
occupy unique niches". Clearly, entomopathogenic fungi and 
their mycoparasites would fall into the unique-niche category. 
Currently, arthropods are considered to be the most diverse 
and species-abundant group of organisms on the planet with 
estimates of 5—10 million species (Zdegaard 2000). If each 
arthropod species hosts atleast one unique fungal pathogen — 
as postulated for the beetle-infecting Laboulbeniales (Bass 
& Richards 2011) — then entomopathogenic fungi and their 
mycoparasites would constitute an immensely richer group 
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than even the most recent data suggest. Certainly, preliminary 
evidence from studies of the zombie-ant fungi — in which spe- 
cies complexes have been identified that may be composed 
of hundreds of taxa — supports this assumption (Araüjo et al. 
2018, 2020). 


Availability of data and material 


The annotated genomes are deposited in GenBank: accession 
numbers JADHZA000000000 (Torrubiellomyces zombiae) and 
JAFEME000000000 (Niveomyces coronatus). The genome as- 
sembly, gene predictions and functional annotations can also 
be analysed interactively at https://fungalgenomics.science. 
uu.nl. Sequences generated in this study for phylogenetic 
analysis have also been deposited in GenBank (see Table 1 
for accession numbers). Holotypes are deposited at the New 
York Botanical Garden Herbarium (type numbers NY4434800 
and NY4434801 for N. coronatus and T. zombiae, respectively). 
Cultures are deposited in the culture collection of the Westerdijk 
Fungal Biodiversity Institute (CBS 149186 and CBS 149187 for 
N. coronatus and T. zombiae, respectively). 
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